Hi-GAL: the Herschel infrared Galactic Plane Survey 

S. Molinarii'2, B. Swinyard^^ J. Bally^ M. Barlow^ J.-P. Bernard^ P. Martin^^, T. 
Moore'^^ A. Noriega-Crespo^^ R. Plume'^^ L. Testi2^'^^ A. Zavagno^ 

and 

o 

O A. AbergeP, B. Ali^, P. Andre^ J.-P. Baluteau^, M. Bcncdcttini^ O. Berne^°, N. P. 

^ Bmot^\ J. Blommaert^^ S. Bontemps^'^^ F. Boulanger^^ J. Brand^l C. Brunt^^, M. 

^ Burton^'', L. Campeggio^^, S. Carey^^, P. Caselli^°, R. Cesaroni^^, J. Cernicharo^^, S. 

I— * Chakrabarti^^, A. Chrysostomou^^'^^ C. Codella^^ M. Cohen^^, M. Compiegne^'^, C. J. 

m Davis^^, P. de Bernardis^^, G. de Gasperis^^, J. Di Francesco^", A. M. di Giorgio^, D. 

^ Elia^^ F. Faustini^, J. F. Fischera^^ Y. Fukui^^ G. Fuller^^^ K. Ganga'^^ P. Garcia-Lario^^ 
r-n M. Giard^ G. Giardino^^, J: Glenn^, P. Goldsmith^^, M. Griffin^s, M. Hoare^o, M. 

^ Huang39, B. Jiang^o, C. Joblin^, J. Joncas^S M. Juvela^^ J. Kirk^s, G. Lagache^^^ J. Z. 

P Li^^ T. L. Lim^^ S. D. Lord^, P. W. Lucases, B. Maiolo^^ M. Marengo^^, D. Marshall^ S. 
^ Masi^^, F. Massi^^ M. Matsuura^'^^, C. Meny^, V. Minier^ M.-A. Miville-Deschenes^, L. 

^ Monticr^, F. Mottc^, T. G. Miiller^^ P. Natoli^^, J. Neves^^, L. Olmi^i, R. Paladini^^ D. 

O Paradis^°, M. Pcstalozzi^, S. Pezzuto^, F. Piacentini^^, M. Pomares'^, C. C. Popescu"^^, W. 
^ T. Reach^ J. Richer^^, I. Ristorcelli^ A. Roy^^ P. Royer^^ D. Russeil^ P. Saraceno^, M. 

Sauvage^ P. Schilke^\ N. Schneider-Bontemps^ F. Schuller^\ B. Schultz^ D. S. 
' ' Shepherd^^ B. Sibthorpc^^ H. A. Smith^^ M. D. Smith^^ L. Spinoglio^ D. Stamatellos^^ 
^ F. Strafellal^ G. Stringfellow^ E. Sturm^^ R. Taylor5^ M. A. Thompson's, R. J. TuffsS^ 
^ G. Umana^^, L. Valenziano^^, R. Vavrek^^, S. Viti®, C. Waelkens^', D. Ward-Thompson^^, 
O G. White^^'5^ F. WyrowskiS^ H. W. Yorke^^, Q. Zhang^o 

(N 
O 

O ABSTRACT 

'. 1 Hi-GAL, the Herschel infrared Galactic Plane Survey, is an Open Time Key Project of the 

^ Herschel Space Observatory. It will make an unbiased photometric survey of the inner Galactic 

k> Plane by mapping a two-degree wide strip in the longitude range | / |< 60° in five wavebands 

between 70/xm and 500/Ltm. The aim of Hi-GAL is to detect the earliest phases of the formation of 

^ molecular clouds and high-mass stars and to use the optimum combination of Herschel wavelength 

coverage, sensitivity, mapping strategy and speed to deliver a homogeneous census of star-forming 
regions and cold structures in the interstellar medium. The resulting representative samples will 
yield the variation of source temperature, luminosity, mass and age in a wide range of Galactic 
environments at all scales from massive YSOs in protoclusters to entire spiral arms, providing an 
evolutionary sequence for the formation of intermediate and high-mass stars. This information is 
essential to the formulation of a predictive global model of the role of environment and feedback 
in regulating the star-formation process. Such a model is vital to understanding star formation 
on galactic scales and in the early Universe. Hi-GAL will also provide a science legacy for decades 
to come with incalculable potential for systematic and serendipitous science in a wide range of 
astronomical fields, enabling the optimum use of future major facilities such as JWST and ALMA. 

Subject headings: ISM-star formation-high-mass stars-IR-Herschel 
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1. Introduction 

Dust is the most robust tracer of the 'Galactic 
ecology' - the cycling of material from dying stars 
to the ionized, atomic, and molecular phases of 
the ISM, into star forming cloud cores, and back 
into stars. While atoms, ions, and molecules are 
imperfect tracers because they undergo complex 
phase changes, chemical processing, depletion onto 
grains, and are subject to complex excitation con- 
ditions, dust is relatively stable in most phases of 
the ISM. It is optically thin in the Far Infrared 
(FIR) over most of the Galaxy, so that its emis- 
sion and absorption simply depend on emissivity, 
column density and temperature. Gold dust in 
particular (10K<T< 40 K) traces the bulk of non- 
stellar baryonic mass in all of the above "habitats" 
of the Galactic ecosystem. 

Temperature and luminosity and, as their by- 
product, mass of cold dust measured over the en- 
tire Galactic Plane (GP), are, at sub-parsec reso- 
lution, the critical quantities needed to formulate 
a global predictive model of the cycling process be- 
tween the Galactic ISM and star formation. This 
process drives the galactic ecology in normal spi- 
rals as well as the enhanced star-formation rates of 
starburst galaxies and mergers and a quantitative 
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understanding of it is needed in order to follow 
the formation and evolution of galaxies through- 
out the cosmos. The adequate measurement of 
these key quantities has been beyond the capabili- 
ties of the previous mid- to far-infrared surveys of 



the Galactic Plane (IRAS, Neugebauer et al. 1984 
MSX, [ Price et al.||200l| COB E/DI RBE and FI 



RAS, e.g. ISodroski et al.|[T994l ISO, |Omont"erar 



2003 



Spitzer, Benjamin et al.|[2003[ [Carey et al. 



2009) either due to limited wavelength coverage 
and/or inadequate spatial resolution leading to 
confusion. The balloon-borne BLAST experiment 
( [Pascale et al]|2008[ ) implements Herschel/SPIRE 
detector arrays and is providing exciting anticipa- 
tions of what Herschel will do. The AKARI satel- 



lite (Murakami et al. 2007) improves over IRAS, 



and results from its FIR photometric mapping of 
the GP are eagerly awaited. 

Observing the distribution and temperature of 
dust across the Galaxy will resolve many current 
debates such as the modes of formation of molec- 
ular clouds and high-mass stars. 

Molecular clouds are traditionally thought to 
follow a "slow formation" scheme, where dis- 
tributed material is accumulated by large-scale 
perturbations such as the passage of a spiral arm. 
Shielding by dust and surface reactions on grains 
promotes the Hi— >-H2 transition, which in turn 
allows the formation of other molecules that cool 
the cloud. Gravity, mediated by magnetic fields, 
leads to star formation. In this scenario cloud life- 
times are about ~30 Myr ( Leisawitz et al.|[l989 |. 
This picture has difficulty explaining the absence 
of quiescent, non star forming GMCs (however, 
Palla & Galli 1997) and the continuous re- 



see 



generation of turbulence needed to support GMCs 
for many crossing times. Alternatively, a "fast for- 
mation" scenario has been proposed (Hartmann 



et al. 2001) in which most MCs are transient. 



short-lived structures (Stone et al. 1998 Padoan 



& Nordlund 1999) created in the post-shock re- 



gions of converging large-scale flows. Stars form on 
very short timescales ( Elmegreen|2000 ) . However, 
rapid MC formation requires rapid Hi— 7>H2 con- 
version (Goldsmith fc Li 2005). Accelerated H2 
formation requires either high-density pre-shock 
conditions (n~200 cm'^, T<100K; [Price et al 



2001), or strong turbulence (Glover & Mac Low 



2007), higher than observed. 

On the other hand, the formation of high-mass 



stars and of the star clusters hosting them is likely 
the most important process that shapes the forma- 
tion and evolution of galaxies. Massive stars are 
responsible for the global ionization of the ISM. 
Their energetic stellar winds and supernova blast 
waves direct the dynamical evolution of the ISM, 
shaping its morphology, energetics and chemistry, 
and influencing the formation of subsequent gen- 
erations of stars and planetary systems. Despite 
their importance, remarkably little is known about 
how massive stars form ( [McKee fc Tan|2003[ ). We 
lack a "fundamental theory" or, rather, a galaxy- 
scale predictive model for star formation. One of 
the main limitations to this goal is the lack of sta- 
tistically significant and well-characterized sam- 
ples of young massive stars in the various evolu- 
tionary stages and environments on which a theory 
can be based. In turn, this results from the diffi- 
culty of gathering observational data on on a large 
number of forming high-mass stars: they make up 
only a very small fraction of the total number of 
stars in the Galaxy, their early evolutionary phases 
of massive stars are more rapid than those of low- 
mass stars, they lie at large distance and form in 
crowded environments. 

There is thus a long list of questions that the 
community has been addressing for some time, not 
finding satisfactory answers. Here is an abridged 
list: 

• What is the temperature and density structure 
of the ISM? How do molecular clouds form, 
evolve, and how are they disrupted? 

• What is the origin of the stellar initial mass 
function (IMF)? What is its relationship to the 
mass function (MF) of ISM structures and cloud 
cores on all scales? 

• How do massive stars and clusters form and 
how do they evolve? What are the earliest 
stages of massive star formation and what are 
the timescales of these early phases? 

• How do the Star Formation Rate (SFR) and Ef- 
ficiency (SFE) vary as a function of Galacto- 
centric distance and environmental conditions 
such as the intensity of the Interstellar Radi- 
ation Field (ISRF), ISM metallicity, proximity 
to spiral arms or the molecular ring, external 
triggers, and total pressure? 
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Fig. 1. — Left panel: / — b plot of A-rising SED IRAS sources; the straight thin line is the Galactic midplane. 
The asterisks mark the median latitude of the sources computed in 10°l bins. The thin lines delimit the 
regions where 50% and 75% of |6| < 5° IRAS sources are contained. Right panel: 5-distribution of the same 
IRAS sources in the |^| < 60° region. 



• Does a threshold column density for star forma- 
tion exist in our Galaxy? What determines the 
value of this possible threshold? 

• What are the physical processes involved in trig- 
gered star formation on all scales and how does 
triggered star formation differ from spontaneous 
star formation? 

• How do the local properties of the ISM and the 
rates of spontaneous or triggered star formation 
relate to the global scaling laws observed in ex- 
ternal galaxies ? 

Using the Herschel telescope, the largest ever 
in space, Hi-GAL, the Herschel infrared Galactic 
Plane survey, will provide unique new data with 
which to address these questions. Hi-GAL will 
make thermal infrared maps of the Galactic Plane 
at a spatial resolution 30 times better than IRAS 
and 100 times better than DIRBE, from which a 
complete census of compact source luminosities, 
masses, and spectral energy distributions (SEDs) 
will be derived. Source distances are a crucial pa- 
rameter in this respect, and a dedicated effort will 
be needed (see Q. Extraction of statistically sig- 
nificant samples of star-forming regions and cold 
ISM structures will be possible in all the environ- 
ments of the Milky Way at all scales from massive 
Young Stellar Objects (YSOs) in individual pro- 
toclusters to complete spiral arms. 



In the following we present the specific charac- 
teristics of the survey as well as some of the sci- 
ence outcomes that we expect to obtain with this 
unique project. 

2. Hi-GAL Observing Strategy 

The area covered by Hi-GAL (| / |< 60°, | 6 |< 
1°) contains most of the star formation in the 
Galaxy, and it is the one which offers the best 
coverage in ancillary data which will be critical in 
the scientific analysis (see The b distribution 
and extent of the survey is shown in Fig. [T] along 
with the I — b plot of A-rising SED IRAS sources 
(-F'loo > Pqo > F25 > F12) which are potential 
YSOs. The Hi-GAL area (thick dashed lines in 
that figure) represents the |&| <1° strip centered 
on the midplane and contains about 80% of all 
potential YSOs contained in |5| < 5° strip, thus 
encompassing most of the potential star formation 
sites in the inner Galaxy. 

The Herschel photometric cameras PACS 



et al. 2008) and SPIRE (GrifHn et al. 2009 



Pogli ;sch 



will 



be used in parallel mode (pMod^ to maximize 
survey speed and wavelength coverage. Due to the 
instruments wavelength multiplexing capabilities. 



^In pMode the Herschel telescope is scanning the sky in 
a raster fashion at constant speed while both PACS and 
SPIRE acquire data simultaneously 
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one pMode observation delivers maps at five dif- 
ferent wavelengths: 70 and 170/im with PACS and 
250, 350 and 500/im with SPIRE. Both cameras 
cameras use bolometric detector arrays to map the 
sky by scanning the spacecraft along approximate 
great circles. Both instruments require their on 
board sub-kelvin coolers to be recycled to provide 
the detectors with the operating temperature re- 
quired of about 0.3 K in each case. In pMode both 
instruments are placed into their photometric ob- 
serving mode with the detectors at their correct 
operating temperature, i.e. both instrument cool- 
ers are recycled, and data are taken from the five 
arrays simultaneously as the spacecraft is scanned 
across the sky. 




Fig. 2.— The field of view of the PACS and SPIRE 
instruments shown in the context of the Herschel 
field of view as viewed in the co-ordinate system 
of the spacecraft +Z refers to the axis towards the 
Sun. The +X axis is along the telescope bore- 
sight out of the page. The scan directions used to 
map the sky are as indicated. The different photo- 
metric channels of each instrument map the same 
region of the sky. 

The size and separation of the fields of view of 
PACS and SPIRE are shown in Fig. [2] as viewed 
in the spacecraft co-ordinate system. Although 
the PACS array fully spatially samples the point 
spread function from the telescope it still has gaps 



between the sub-arrays, and the SPIRE arrays 
only sparsely sample the sky. In order to make 
fully spatially sampled maps it is necessary to scan 
the SPIRE array at an angle of 42.5° with respect 
to its short symmetry axis. Scanning at an angle is 
also used for the PACS arrays to fill in for the gaps 
between sub- arrays. To achieve redundancy in the 
data and remove instrumental effects such as high- 
frequency detector response, slow drifts in gain or 
stray light, saturation and environmental (cirrus 
confusion) effects it is also necessary to make at 
least a second pass over the same region of the sky 
using the other scan angle at —42.5° angle which, 
quite conveniently, is nearly orthogonal to the first 
one (see fig. [3|. 

The distance between each scan in parallel 
mode is set by the size of the PACS array (being 
the smaller of the two), and the effective length of 
each leg of the raster takes into account the sep- 
aration between the two fields of view. HSPOT, 
the Herschel-SPOT observing tooQ automatically 
calculates these parameters to ensure that the area 
required is covered. The distance between scans is 
approximately 155" and the excess length of the 
scan beyond the required length to cover the area 
is typically 20'. An example of how the sky is cov- 
ered in a Parallel Mode observation used in Hi- 
GAL is shown in Fig. [3) 




Fig. 3. — Sample AORs (Astronomical Observa- 
tion Request) overlaid on the IRAS lOO/xm im- 
age of a portion of the Galactic plane. From left 
to right we show: both nominal and orthogonal 
2.2°x2.2° pMode AORs overlaid on one another. 
The PACS-covered area is outlined in pink while 
SPIRE is green. 



^ftp;/ /ftp. sciops.esa.int/pub/hspot/HSpot_download. html 
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The strategy employed to cover the — 60°< I < 
60°, |6| <1° survey area is to conduct observations 
with series of 55 2.2° x2.2° square tiles with two 
passes over each tile with the two above mentioned 
scan angles (see Fig. [2]). These tiles will be spaced 
every 2° , so that the overlap between tiles ensures 
that no coverage gaps are introduced by different 
tile orientation due to variable satellite roll angles 
with time; fig. |3] shows a section of the galactic 
plane with consecutive observing blocks overlaid 
providing overlapping coverage. 

Given the spatial separation required for PACS 
in the pMode observations, the SPIRE data is 
heavily oversampled and we cover a greater area 
than required with each individual instrument 
than would be required using them sequentially. 
Although it might seem that sequential PACS and 
SPIRE scan mode observations would be more ef- 
ficient in fact the satellite overheads, set up, cal- 
ibration and pointing acquisition, et cetera, mean 
that it requires 30% more time to cover the same 
area sequentially compared to using the pMode. 

In order to cover the maximum area in the 
shortest time Hi-GAL data will be taken at the 
maximum possible scan speed for the satellite of 
60 "/sec. This implies a beam crossing time for 
the short wavelength, 250 /im band of SPIRE of 
3 Hz well within the bandwidth available in the 
detectors of 5 Hz. However, although the PACS 
detectors have a similar response time the much 
smaller point spread function will be smeared out 
compared to that achievable with a slower scan. 
Additionally because of the finite data transmis- 
sion bandwidth between the Herschel satellite and 
the ground, it is necessary to perform on-board 
data compression for the PACS data which are 
the most demanding in terms of number of pixels 
(2048 for the 70/im array and 512 for the 170/im 
array) at the frame acquisition rate of 40Hz. The 
baseline configuration for the pMode is then to av- 
erage on board groups of 8 frames at lOfim and 4 
frames at 170/ini. Since the telescope is contin- 
uously scanning while acquiring, this coaddition 
will result in a further degradation of the Point 
Spread Function in the direction of the scan from 
its original diffraction limited shape; the effect will 
be more severe at 70/im where the degradation 
should be of a factor two based on simulations. 
This loss in imaging fidelity at the shortest wave- 
length is considered acceptable for a survey like 



Hi-GAL because, as discussed in the previous sec- 
tion, our main focus is toward a large-scale picture 
of the galaxy. Taking advantage of the orthogonal 
cross scan observing strategy, we may be able to 
recover some of the spatial resolution by careful 
deconvolution during post processing. 

2.1. Detection of compact sources 

The SPIRE digital readout electronics impose 
a limitation on the brightest sources that can be 
observed for a given offset setting (DC voltage 
removal) before digitization (SPIRE Instrument 
Users Manual, 2007). This problem can be al- 
leviated to some extent by choosing a bias set- 
ting that gives the largest dynamic range per offset 
range. Simulations of the effect of bias variation 
show that setting a bias higher (~3x) than the 
predicted nominal value will approximately dou- 
ble the dynamic range for most detectors under 
the conditions likely to be found in orbit (tele- 
scope temperature and emissivity and sky back- 
ground). The same simulations show that a sig- 
nificant (>10%) fraction of the SPIRE 250^m ar- 
ray detectors will saturate on sources greater than 
500 Jy. The situation is slightly more relaxed for 
the 350 and 500/im arrays. We take the upper 
limit of detectable sources in the SPIRE bands as 
500 Jy/beam assuming that a strong source in- 
strument setting is used. This setting is required 
for all observations of bright regions/sources with 
SPIRE and is not a special Hi-GAL configuration. 
The saturation limits for PACS should be around 
2000Jy at nominal bias, that will be used for the 
Hi-GAL survey. 

The l-fj sensitivities provided by HSPOT, for a 
single Astronomical Observation Request (AOR), 
are 17.6 and 26.8 mJy in the two PACS 70 and 
170/im bands, and 12.8, 17.6 and 14.9 mJy for 
the SPIRE bands; co-addition of the orthogonal 
scanning patterns will provide \/2 better figures. 
These sensitivities result from the adopted scan- 
ning strategy designed to maximize redundancy 
and map fidelity especially for large scale diffuse 
structures. However, the limiting factor for the de- 
tectability of sources and clouds will likely be cir- 
rus confusion. Estimates based on recent BLAST 
measurements ( Roy et aL]|2010 l suggest values of 
the order of 75, 140, and 160 mJy in the 170, 250 
and 500/im Herschel bands for a representative re- 
gion of the Galactic Plane at ^=45°; these values 
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Fig. 4. — Flux expected in the 170, 250 and 
500/xm Herschel bands (color coded) from a 2OM0 
core (dust+gas) for different temperatures (differ- 
ent symbols), as a function of distance in kpc. 
The dashed lines represent the confusion noise ex- 
pected in the three bands. 

are greater than the detector sensitivities. Fig. |4] 
shows the expected flux from a 2OM0 envelope as 
a function of distance (in kpc) in each of three 
bands above mentioned, and for three different 
dust temperatures; we adopted (3=2 and the dust 
opacity from Preibisch et al. (1993). The hori- 



zontal dashed lines (color-coded with wavelength) 
represent the predicted confusion noises based on 
BLAST images. Fig. |4] shows that we will de- 
tect the representative 2OM0 core everywhere in 
the Galaxy except for very cold dust (T<10K), 
for which detectability is predicted to be limited 
within a distance of about 5 kpc. We may then 
conclude that cirrus confusion is not going to be 
a problem for the investigations of the interme- 
diate and high-mass star formation studies which 



are the "core science" of this project (see ^3.3 and 
311). 



2.2. Detection of extended structures 

The diffraction-limited instrument beams at all 
wavelengths can be used to translate the confusion 
noises reported in the previous section into bright- 
ness units to investigate the detectability limits 
expected for diffuse extended structures. Calcula- 
tions of the expected brightness levels from opti- 
cally thin dust as a function of temperature and H 
column density (assuming gas/dust=100) are re- 
ported in Fig. [5^,b for A=170 and 500/xm respec- 




10^' 10^- 

N[H] (cm" 



Fig. 5. — Brightness (in MJy/sr) of optically thin 
dust as a function of temperature and column den- 
sity at 170/im (top) and 500/im (bottom). The red 
lines indicate the expected confusion noise. 



tively, where the expected confusion noise levels 
are also reported in red lines. 

The figures show that with a typical confusion 
noise of about 10 MJy/sr at 500/im it will be pos- 
sible to detect most Infrared Dark Clouds (IRDCs, 
see ^3.2), where column densities are in the range 
lO^'^Hr^cm"^. Less dense clouds with column 
densities of the order of few lO^^cm"^ should be 
easily detectable at levels of tens of MJy/sr at 
170/im at temperatures as low as T~20K. The 
situation may be less simple in the regions closer 
to the Galactic center. However, our broad spec- 
tral coverage provides an important advantage 
for measuring the temperature accurately, and 
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for isolating structures and sources with tempera- 
ture different from the standard diffuse ISM cir- 
rus (^20K). Besides, all of the above is based 
on extrapolated estimates and we will provide the 
definitive measurement of the cirrus confusion at 
Herschel wavelengths and resolution. 

2.3. Checks on photometric calibration 

Virtually every block of the Hi-GAL survey will 
contain secondary calibrators, either stars or other 
well known objects, ensuring accurate checks of 
the flux calibration of Hi-GAL data. Our base- 
line calibration targets will be the 400 stars used 
as calibrators for the Spitzer GLIMPSE-I/II sur- 
veys. These are mostly A0-5V or KO-MOIII stars, 
although they also include ~60 calibrators of hot 
or warm dwarfs (B-G). In order to remove sources 
with non-photospheric FIR emission we will make 
predictions for the 24/Ltm fluxes and then test for 
excesses in the MIPSGAL data first. Once anoma- 
lous sources are excluded we can extrapolate to the 
FIR and create an initial set of calibrators. 

We will also be able to obtain a reliable calibra- 
tion for extended sources, which is one of the more 
difficult parts of the nominal instrumental calibra- 
tion activity. As part of the calibration scheme 
we will compare fluxes in the SPIRE 500/xm band 
with fluxes from the same band of Planck-HFI in 
suitable locationtlf] 

3. Hi-GAL Key Science 

3.1. The Distribution of the ISM Temper- 
ature and the Intensity of the Inter- 
stellar Radiation Field 

At near-infrared wavelengths, the emission 
from dust is produced by small particles whose 
abundance varies significantly, being strongly de- 
pleted by coagulation processes in the dense ISM 
(e.g. [Bernard et al.|[l993| [Abergel et"aI1|1994 
Stepnik et al.||2003[ ). Far-infrared (FIR) emission 



is produced by larger grains which are more sta- 
ble and dominate the total dust mass and trace 
all phases of the ISM. The ISM dust spectrum 
peaks in the FIR where the Galaxy is transpar- 
ent. FIR emission is therefore a reliable tracer 
of the overall ISM column density structure in 

•^many of the Hi-GAL Co-Is are also Planck Consortium 
members 



our Galaxy. Other phase-independent tracers in- 
clude dust absorption and gamma-ray production, 
where however the former can be used to sample 
only the nearest 1 kpc, and 7-ray surveys currently 
lack sensitivity and angular resolution. 

Variations in the FIR emissivity (the ratio of 
surface brightness to column density) are domi- 
nated by the non-linear effects of dust temper- 
ature through the Planck function. Fortunately, 
the shape of the dust SED as measured by PACS 
and SPIRE will be most sensitive to temperature 
variations as the spectral bands sample the peak 
of the Big Grain emission and the contribution 
of Very Small Grains can be estimated from the 
Hi-GAL data at 70^m and MIPSGAL at 24^m. 
The dust temperature (Td) and its spatial varia- 
tions will therefore be measured precisely. This 
important parameter can be used, in conjunction 
with complementary data from Planck, HI, CO, 
Ha and 7-ray surveys , to estimate the strength 
and spectral shape of the InterStellar Radiation 
Field (ISRF), which is set by the stellar content 
in a given region. So far, the dust temperature in 
the Galactic Plane has been mapped over limited 



regions using IRAS (Kim et al. 1999 Douglas & 



|Taylor|2007"| and over the entire Plane at a resolu 
tion of 40' ( |Lagache et aT[l998) using DIRBE. Hi- 
GAL will improve with respect to the latter by a 
factor of about 100 in linear scales. It will trace the 
local radiation field on scales relevant to star for- 
mation, and provide mass estimates even at large 
distances. In the case of the dense medium, deter- 
mining the 3-D distribution of the ISRF strength 
and spectral shape in a given cloud will require 
radiative transfer modeling. This is possible, even 
for complex geometries, using Monte-Carlo codes 



(e.g. Juvela & Padoan 2003). Using such codes. 



the equilibrium dust temperature Td and the dust 
emission can be predicted at any 3D location in 
the cloud. Integration along the line-of-sight in 
turn allows to predict 2D emission maps. 

3.2. Molecular cloud formation 

About a quarter of the mass in the ISM is in 
molecular form ( Blitz|[l997 ) and most of that ma- 
terial resides in giant molecular clouds (GMCs). 
Since GMCs are also the dominant sites of star for- 
mation, understanding their origins and evolution 
is essential to our understanding of the Galactic 
environment. 
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Fig. 6. — Model of a typical IRDC after Stamatellos et al. (2004) at seven different wavelengths illustrates 



the importance of the Herschel wavelength coverage and resolution (at 2kpc the linear extent is 5') to obtain 
measurements of IRDCs in emission. 



In combination with molecular line surveys, Hi- 
GAL will provide the data needed to derive basic 
physical properties of GMCs. We will detect and 
characterize cold structures in the inner GP and 
classify them based on star formation activity. De- 
tection statistics for clouds with different temper- 
atures and degree of star formation activity will 
provide the fraction of quiescent vs star-forming 
clouds. It will thus be possible to constrain the 
properties and lifetimes of GMCs in our Galaxy 
and to compare with the predictions of fast evo- 
lution of molecular clouds ( Hartmann et al.||2001 ) 
or a more traditional slow evolution of star for- 
mation in our Galaxy (Shu et al. 1987). Varia- 



of insufficient wavelength coverage (e.g. IRAS) or 
inadequate spatial resolution (DIRBE, FIRAS). 
CO observations are problematic due to molecu- 
lar freeze-out onto grains ( [Flower et aT1[2005| , or 
photo-chemical effects in low-metallicity environ- 
ments (Bot et al. 2007). The recent detection of 



very cold clumps in the GP with Archeops ( Desert 



tions with Galactocentric radius will determine if 
the slow/fast scenarios are mutually exclusive or 
reflect different initial/environmental conditions. 
A large-area survey like Hi-GAL will provide the 
needed statistical significance in all mass bins, es- 
pecially at the high-mass end, and in a variety of 
Galactic environments. 

Direct detection of cold (i.e. T<20 K) dust 
which could be the quiescent counterparts to 



et al. 20081, confirms the FIR and submm con- 
tinuum as the best tool to trace cold ISM compo- 
nents. Notable examples are Infrared Dark Clouds 
(IRDCs) and Hi Self- Absorption (HISA) clouds. 

IRDCs are structures initially discovered as ex- 
tinction features against the bright mid-IR Galac- 
tic background, and soon verified to exhibit prop- 
erties similar to molecular clouds. Their prop- 
erties (n > lO^cm-^^ Nh ^ lO^^ - lO^" cm^^ 
and T < 25K - | Egan et al.] [19981 [Carey et al. 



2000 



2006 



i? ~ 5pc and M ~ lO'^ Me Simon et al 



Rathborne et al.|2006 1 suggest that they are 



GMCs, is diflicult ( 


Sodroski et al. 


1994 


Reach 


et al. 


1995 


Lagache et al. 


1998 


1 either because 



the precursors of cluster-forming molecular clumps 
like Orion. Thus, IRDCs are ideal for the study of 
the pristine, undisturbed physical conditions that 
may produce massive stars and clusters. IRDCs 
have only been detected against the bright Galac- 
tic mid-IR background (mostly for |^| < 30°); their 
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true Galactic distribution is unknown. Modeling 
of IRDCs in the IR and submm (Fig. [6| [Stamatel 
los et al.|2004 ) proves Herschel's unique ability to 
detect them and to measure their SEDs. IRDCs, 
with T200fj.m > 1, are not detectable by either 
IRAS or Spitzer. Hi-GAL will provide a defini- 
tive inventory of cold dust and potential sites of 
massive star cluster formation everywhere in the 
inner GP. 

HISAs are traced by cold Hi gas seen in absorp- 
tion against a background of warm Hi emission 
dGibson et aL][2000l [20051 [Goldsmith fc"Li||2005J ; 
they may provide additional clues to the formation 
of molecular clouds. When compared with molec- 
ular tracers, these cold (relative to the ambient 
neutral medium) and relatively quiescent {Av ~ 
1 — 3kms) clouds show a wide range of H1/H2 ra- 



tios (Li & Goldsmith 2003 Klaassen et al. 2005) 



which suggests that they might be Hi— >-H2 con- 
version sites. Detailed studies of their FIR— ?>mm 
SED shapes can help clarify this issue, providing 
evidence for the grain types necessary for a reason- 



able Hi— ^H2 formation timescale ( Goldsmith & Li 



2005). 



3.3. Timeline of high-mass star formation 

The paradigm for the formation of solar-type 
stars via accretion through a circumstellar disk 
( Shu et ar][1987 ) predicts an evolution from cores 
to protostars and, finally, pre-main sequence stars 
that is well matched with distinctive characteris- 



tics of their SEDs (Lada & Wilking 1984 Andre 



et al. 1993). The empirical classification of the 
SED of low mass YSOs has thus been used as a 
powerful tool to constrain theoretical models. 

Higher mass stars reach the conditions for H- 
burning faster than the time required to assemble 
them, so that winds and radiative feedback will 
strongly influence accretion and may limit the fi- 



nal mass of the star (Zinnecker & Yorke 20071. 
However, since massive stars exist, several theories 
have been proposed to solve this puzzle including 
accretion from a disk, a very high-pressure ambi- 



ent medium, "flashlight effect" 


Yorke & Sonnhal- 


ter 


2002 


McKee & Tan 


2003 


, competitive ac- 


cretion, or coalescence (e.g. 


Stabler et al. 


2000J. 



Application of SED-based classification tools, and 
evolutionary diagnostics like the Menv-^boi dia- 
gram which relates the bolometric luminosity of 



2008 ), to a large sample of luminous protostar can- 
didates in the GP will define a timeline for the 
various phases of massive star formation that will 
constrain the theories and lead to new estimates of 
the SFR. Clearly the source distance is the crucial 
parameter here; we are collecting the information 
from the major molecular line surveys over the in- 
ner Galactic Plane, while planning to undertake 
additional surveying activities at a variety of fa- 
cilities in several high-density tracers to get addi- 
tional data (see J|4[for more detail). 

An evolutionary sequence for massive YSOs 
has been proposed (cold massive cloud core; Hot 
Molecular Core with outflow; IR-bright massive 
YSO; ultracompact (UC) Hii region, e.g. ,Evaiis[ 
et al.[[2002[ [Kurtz et al.][2000[ ) but it is qualitative 
and based on small and possibly incomplete sam- 
ples. Samples of bright and massive YSOs (Moli- 



nari et al.|1996[|Sridharan et al.|2002l|Hoare et al 



2004[ ) are IRAS or MSX selected and tend to suf- 
fer from age biases and confusion which prevent 
firm quantitative conclusions. A phase of intense 
and accelerating accretion prior to H-burning igni- 



tion, that may be observable (e.g. Molinari et al 



1998) in the form of dense condensations devoid 



of IR as well as radio continuum emission, seems 



confirmed by recent large mm surveys (Beltran 



et al. 


2006 


Hill et al. 


2005 



Millimeter contin- 
uum alone, however, cannot distinguish between 
pre-collapse condensations and rapidly accreting 
cores; Hi-GAL will use the full potential of Her- 
schel wavelength coverage and spatial resolution to 
trace the SED peak of dust envelopes in all phases, 
from massive pre-stellar condensations to UCHii 
regions. An angular resolution of 30" or less, 
typical of Hi-GAL, has been proven (e.g 



Moli- 



a YSO to the mass of its envelope (Molinari et al 



nari et al. 2008) to be the key to building accu- 
rate SEDs, deriving reliable luminosities for mas- 
sive YSOs, and distinguishing embedded UCHiis 
strongly emitting in the Mid-IR and radio from 
pre-UCHiis (Fig. ^. 

The abundance of high-mass YSOs per mass 
bin in the various evolutionary phases will provide 
an estimate of the duration of each phase. This 
timeline can be directly compared with the pre- 
dictions of various models, and together with the 
YSO mass function, will be used to infer the SFR. 
As an example, using current estimates for SFR 
and IMF jMcKee fc Williams|1997[[Kroupa|200l"| ), 
and a 10^ yr period to assemble a massive star 
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60" 




Fig. 7.— Mid-IR 21^m MSX images for a candidate precursor of a Hot Core (IRAS23385+6053, left), 
and a UCHii (IRAS05137+3919, right). Full and dashed contours represent the millimeter (cold dust) and 

Hi- GAL will complement these with similar 



radio (jet) continuum emission (from Molinari et al 
resolution FIR data. 



(2008)) 



( |McKee fc Tan|2003[ ) , we expect > 1000 FIR/sub- 
mm objects with M>15M0 in the Galaxy, and 
only --180 objects with M>50Mo. Recent mil- 



limeter surveys in the Cyg X region ( Motte et al 



2007 ) confirm the rarity of such massive precursors 



and strengthen the need for a systematic unbiased 
search. 

3.4. Bridging the Gap between Global and 
local star formation 

Galactic phenomenology currently invokes an 
indeterminate mixture of spontaneous and trig- 
gered star formation. Triggering agents include 



radiation pressure from OB stars ( |Sugitani et al 
1989), compression by expanding Hll regions 
( Elmegreen fc Lada|[l977l [Deharveng et al.|[2005| , 
or fragmentation of supershells by multiple su- 



pernovae in OB associations (McCray & Kafatos 



1987). On larger scales we still do not know if 



spiral density waves actively induce star forma- 
tion (e.g., Elmegreen 2002) or simply assemble 



star-forming regions, with local feedback and trig- 
gering becoming more important within the arms 



(Sleath & Alexander 1996). The mean SEE of 



a galaxy can increase (up to 50 times) in star- 



bursts (Sanders et al. 1991) and galaxy mergers 



due to strong feedback effects, a process observed 
in miniature in Galactic star-forming regions (e.g. 
Moore et aL]|2007[ ). Whether the IMF depends on 



local triggering and other environmental factors is 
unclear. 

Hi-GAL will enable quantitative analysis based 
on basic observables - the luminosity functions of 
YSOs, the mass function of dense star-forming 
structures and quiescent clouds. Hi-GAL will pro- 
vide the essential context of high-mass star for- 
mation, as it relates to molecular gas. Hi gas, 
stars, Hii regions, OB associations, SNRs and spi- 
ral arms. Theoretical models and numerical sim- 
ulations will be tested in multiple ways. We will 
discover whether a local triggering agent is nec- 
essary for high-mass star formation or if a spiral 
arm is sufficient, clarifying the differences between 
spontaneous and triggered star formation. We will 
quantify the relationship between the interaction 
strength (estimated using available data from an- 
cillary surveys) and the resulting increase in SEE 



11 



above the spontaneous rate. By locally relating 
the SFR to the properties of the ISM we will probe 
star formation thresholds as a function of envi- 
ronment and spatial scale, and possibly unveil the 
mechanism giving rise to global Schmidt-like scal- 
ing laws. Wc will determine the dominant physical 
process underlying triggering. 

3.5. Serendipitous Science 

The 5-band FIR images and source catalogues 
provided by Hi-GAL will allow research in many 
fields that can only be partially anticipated. A 
detailed description of the specific outcomes in all 
these fields goes beyond the scope of the present 
paper. These aspects will of course be the object of 
dedicated publications to be released in due time. 
In the following we list some of the possible by 
products of Hi-GAL : 

• An input catalogue for ALMA: we expect to de- 
tect some 200-400 objects per tile, most of which 
will mark very cold objects to be studied at all 
possible wavelengths; 

• complete characterization of the Galactic fore- 
ground in the Far-IR and submillimeter, critical 
for the correct interpretation and modeling of 
cosmological backgrounds; 

• Dust formation and destruction in supernovae 
remnants 

• Debris dust disks around main sequence stars, 
with unbiased statistics on frequency and mass 
as a function of star age 

• Evolution of dust properties, especially around 
AGB stars, the factories of cosmic dust; 

• Detection of detached dust shells around first 
ascent giant stars to investigate missing mass in 
AGB envelopes; 

• Detection of multiple shells around AGB stars, 
post-AGB objects and planetary nebulae, as 
well as around various classes of interacting bi- 
naries; 

• Detection of ejecta shells and swept-up ISM 
bubbles around massive stars, providing a com- 
plete census of WR and LBV stars; 

• Extinction maps to aid in correcting Near-IR 
galactic star counts; 



• Detection of Solar system objects via compari- 
son of cross-linked rasters: in particular the de- 
tection of asteroids will be very interesting for 
studies of the albedo 

• Nearby Low-Mass SFRs in the GP: Herschel will 
detect many nearby star forming regions and 
individual YSOs. 

4. Hi-GAL and its place in the context of 
the Multi- Wavelength Milky Way 

PACS and SPIRE are unique in tracing the 
peak of the Spectral Energy Distribution (SED) of 
cold dust and, hence, temperatures and luminosi- 
ties of both star-forming complexes and the ISM, 
at resolutions unmatched by any previous instru- 
ments. However, the full potential of the Hi-GAL 
survey will be realised in the context of the other 
unbiased Galactic Plane surveys which are shaping 
our understanding of the Galactic ecosystem. A 
suite of surveys in the mid- and far-infrared contin- 



uum, ISOGAL (Omont et al. 2003 


I, MSX 


Price 


et al. 


20011, GLIMPSE iBenjamin et al. : 


20031, 


MIPSGAL (Carey et al. 2005), IRAS ( 


Neugebauer 


et al.|1984 


) and Akari (Murakami et al.|2007 


I, has 



been and will be complemented by surveys in the 
submillimeter and millimeter spectral range in- 
cluding the BGPS survey with Bolocam at 1.1mm 
( [Rosolowsky et al.|[2009l ), the ATLASGAL survey 
currently underway with the LaBoCa camera at 



APEX (Schuller et al. 2009), and the SCUBA2 



JPS survey beginning in 2010 (Fig.|8|. 

BLAST and AKARI have wavelength cover- 
age and resolution not too different form PACS 
and SPIRE. However, BLAST has been used to 
map limited portions of the GP while results from 
AKARI photometric imaging of the GP are not 
yet found in the literature. 

Hi-GAL will be the scientific keystone of this 
suite of surveys, completing the continuous cov- 
erage of the dust continuum over three orders of 
magnitude in wavelength at sub-30" resolution, 
and allowing the measurement of dust tempera- 
tures and luminosity over the inner Galactic Plane. 

An extensive plan for radio spectroscopic cross- 
correlation and follow up of the Hi-GAL survey 
has been devised. The top priority is for radio 
spectroscopic observations to obtain distance esti- 
mates for detected sources and structures. CO and 
^■^GO data at sub-arcminute spatial resolution are 
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Fig. 8. — Wavelength-longitude coverage plot for 
photometric (green) and spectroscopic (blue) ex- 
isting or planned GP surveys for the next decade. 
Hi-GAL (red) fills the critical gap where the ISM 
dust emission peaks, between MIPSGAL and the 
SCUBA2/LaBoCa/Bolocam surveys. 

available via the IGPS for the first Galactic Quad- 
rant (QI); lower resolution CO data are available 
from the NANTEN survey (QIV). Should CO res- 
olution not be sufficient, e.g. in case of very cold 
cores with CO depletion, we plan extensive follow- 
up programs (e.g. N2H+, NH3) at the Mopra 22m 
antenna (Australia) for QIV, and at Green Bank, 
Effelsberg, Medicina and Onsala for QI. HI is also 
available for QI and IV at sufficient resolution 
from the IGPS to help resolve the distance am- 
biguity in the inner Galaxy [Busfield et aL||2006l 
Ha VPHAS+IPHAS surveys will be also used. 

In addition to the kinematic distance estimates, 
we plan to use NANTEN2, Mopra and APEX 
for detailed multiline studies of evolution-sensitive 
chemical tracers (N2H+, NH3, CS, HCN, CH3OH, 
CH3CN, etc.) toward clouds and objects discov- 
ered by Hi-GAL. Hi-GAL catalogues will be the 
primary source of major future high spatial reso- 
lution follow ups in the sub-mm with ALMA; in 
the meantime such programs will be attempted us- 
ing the SMA interferometer through Legacy class 
proposals. 



5. Data Processing and Products 

It is relatively easy to translate our scientific 
goals into a clear set of requirements on the data 
processing: we require that the dust continuum 



emission be detectable, and accurately measur- 
able, at all bands over the broadest range in sig- 
nal levels (down to the confusion limit) and spa- 
tial scales. The observing strategy is carefully de- 
signed to that effect, but ensuring that this infor- 
mation is properly extracted from the data stream 
over a 240 sq.deg. area is a formidable challenge. 
We will use the Herschel Interactive Processing 
Environment (HIPE) for all those processing steps 
dealing with fundamental instrumental calibration 
and issues, but we anticipate areas where a ded- 
icated set of specialised tools can take advantage 
of the homogeneous observing strategy and deliver 
higher quality results compared to the standard 
pipeline products; pointing refining, map-making, 
source extraction and photometry are examples. 

We will make available a set of data products 
which will include maps and compact source cat- 
alogues at the five Hi-GAL wavelengths. These 
products will be made available after the end of 
observations (EoO) for the entire survey via in- 
cremental releases. 

Improved reprocessed maps and source cata- 
logues will be subsequently made public, together 
with a first release for an extended source cata- 
logue. All public deliveries will be accompanied 
by an Explanatory Supplement. 

In addition to this minimum set of products, 
quite standard for any large-scale survey like Hi- 
GAL, we plan to make available to the community 
a set of scientific value-added products which will 
be created during our scientific analysis, including 
band-merged catalogues integrated with data from 
continuum surveys at adjacent wavelength, color 
maps, source-subtracted maps. The public access 
to this final set of products is foreseen for EoO-|-42 
months. 

6. Conclusions 

Hi-GAL is an Open Time Key Project to be 
performed with the 3.5 m orbiting Herschel tele- 
scope, to map photometrically the inner Milky 
Way (I I \< 60°,! b \< 1°) in five wavebands 
between 70/xm and 500/im simultaneously, using 
~350 hours of observing time. The unique com- 
bination of survey speed, high sensitivity, high 
spatial resolution and wavelength coverage (across 
the peak of the dust emission) make Hi-GAL the 
first dedicated project to study the early phases 
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of GMC- and high-mass star formation in the 
Galaxy, with a legacy value similar to the IRAS 
mission some 20 years ago. The outcomes of Hi- 
GAL will consist of source lists and images to be 
released in due course after EoO. 

We are grateful to all the people who made 
the building and launch of Herschel such a suc- 
cess. In particular ESA and the Herschel Project 
Scientist G. Pilbratt, and the instrument teams 
of PACS and SPIRE magnificently led by A. 
Poglitsch (MPE, Garching) and M. Griffin (Univ. 
of Cardiff). 

Facilities: Herschel, SPIRE, PACS. 
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